Boosted high harmonics pulse from a double-sided relativistic mirror 
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An ultra-bright high-intensity X- and gamma-radiation source is proposed. A high-density thin 
plasma slab, accelerating in the radiation pressure dominant regime by a co-propagating ultra- 
intense electromagnetic wave, reflects a counter-propagating relativistically strong electromagnetic 
wave, producing strongly time-compressed and intensified radiation due to the double Doppler effect. 
The reflected light contains relativistic harmonics generated at the plasma slab, all upshifted with 
the same factor as the fundamental mode of the incident light. 

PACS numbers: 52.38. Ph, 52.59.Ye, 52.38.-r, 52.35.Mw, 52.27.Ny 



Interaction of electromagnetic (EM) wave with the relativistic mirror has been used by A. Einstein to illustrate basic 
effects of special relativity pQ. In modern theoretical physics the consept of relativistic mirror is used for solving a 
wide range of problems, such as the dynamical Casimir effect [2], the Unruh radiation [3] and other nonlinear vacuum 
phenomena. Relativistic mirrors made by wake waves may lead to an electromagnetic wave intensification [1] resulting 
in an increase of pulse power up to the level when the electric field of the wave reaches the Schwinger limit when 
electron-positron pairs are created from the vacuum and the vacuum refractive index becomes nonlinearly dependent 
on the electromagnetic field strength [5] . In classical electrodynamics EM wave reflected off a moving mirror undergoes 
the frequency and electric field magnitude multiplication, a phenomenon called the double Doppler effect. If the EM 
wave is co-propagating with respect to the mirror, its frequency and energy decreases upon reflection. If it is counter- 
propagating, the reflected light gains energy and becomes frequency-upshiftcd. For the latter case the multiplication 
factor is approximately Aj 2 , where 7 3> 1 is the Lorentz factor of the mirror, making this effect an attractive basis 
for a source of powerful high-frequency EM radiation. Relativistic plasma provides numerous examples of moving 
mirrors which can acquire energy from co-propagating EM waves or transfer energy to reflected counter-propagating 
EM waves (see [5] and references therein) . Manifestation of the Doppler effect in plasma governed by strong collective 
fields is seen in the concepts of the sliding mirror [7J, oscillating mirror [5], flying mirror [3] and other schemes [5] 
which can produce ultra-short pulses of XUV- and X-radiation. 

In this paper we discuss the concept of the accelerating double-sided mirror, Fig. [l] which can efficiently reflect 
the counter-propagating relativistically strong electromagnetic radiation. The role of the mirror is played by a high- 
density plasma slab accelerated by an ultra-intense laser pulse (the driver) in the Radiation Pressure Dominant 
(RPD) regime (synonymous to the Laser Piston regime), | 10| . Such an acceleration can be viewed as the double 
Doppler effect: it is the reflection that allows the energy transfer from the driver pulse to the co-propagating plasma 
slab, which acquires the fraction « 1 - (47) ~ 2 of the driver pulse energy, [10]. The plasma slab also acts as a 
mirror for a counter-propagating relativistically strong electromagnetic radiation (the source). As such it exhibits 
the properties of the sliding and oscillating mirrors, producing relativistic harmonics. The source pulse should be 
sufficiently weaker than the driver, nevertheless it can be relativistically strong. In the spectrum of the reflected 
radiation, the fundamental frequency of the incident radiation and the relativistic harmonics and other high-frequency 
radiation like bremsstrahlung generated at the plasma slab are multiplied by the same factor, w 47 2 , Fig. [T] 
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FIG. 1: (color) The accelerating double-sided mirror concept. 
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Compared with previously discussed schemes, the double-sided mirror concept benefits from a high number of 
reflecting electrons (since the accelerating plasma slab initially has solid density and can be further compressed 
during the interaction) and from the multiplication of the frequency of all the harmonics (since the interaction is 
strongly nonlinear and the mirror is relativistic) . This concept opens the way towards extremely bright sources of 
ultrashort energetic bursts of X- and gamma-ray, which become realizable with present-day technology enabling new 
horizons of laboratory astrophysics, laser-driven nuclear physics, and studying the fundamental sciences, e.g. the 
nonlinear quantum electrodynamics effects. 

In order to investigate the feasibility of this effect we performed two-dimensional (2D) particle-in-cell (PIC) simu- 
lations using the Relativistic ElectroMagnetic Particle-mesh code REMP based on the density decomposition scheme 
[12"] . The driver laser pulse with the wavelength A^ = A = 2nc/uj, the intensity Id — 1.2 x 10 23 W/cm 2 x (1/zm/A) 2 , 
corresponding to the dimcnsionless amplitude = 300, and the duration = 207r /uj is focused with the spot 
size of Dd = 10A onto a hydrogen plasma slab with the thickness I = 0.25A and the initial electron density 
n e = 480n cr = 5.4 x 10 23 cm~ 3 x (1/im/A) 2 placed at x = 10A. The plasma slab transverse size is 28 A. The 
driver pulse shape is Gaussian but without the leading part, starting 5 A from the pulse center along the x-axis. At 
the time t = 0, when the driver pulse hits the plasma slab from the left (x < 10A), the source pulse arrives at another 
side of the slab from the right (x > 10.25A). The driver is p-polarized, i. e., its electric field is directed along the 
y-axis. The source pulse is s-polarized (its electric field is along the z-axis). It has the same wavelength as the driver 
pulse. Its intensity is I s = 1.2 x 10 19 W/cm 2 x (1/zm/A) 2 , corresponding to the dimensionless amplitude a s = 3, its 
duration is r s = 1207r/(jj and its waist size is D s = 20A. The source pulse has rectangular profile along the x- and 
y-axes; such the profile is not necessary for the desired effect but helps to analyse the results. The simulation box 




FIG. 2: (color) (a) The driver and source pulses represented by the y- and ^-components of the electric field, respectively, and 
the ion density (black), (b) The ion energy (curve) and anglular (grayscale) distributions. Both the frames for t — 37 x 2n/u>. 
(c) The electric field 2-component, showing the source pulse overlapped with the first two cycles of the reflected radiation at 
( = 4x 2tv/u. 
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FIG. 3: (color) (a) The electric field component E z representing the reflected radiation along i-axis for t* — 32 x 2n/u>. 
(b) Colorscale: the modulus of the spectrum of E z (r) seen along a;-axis, taken for each t with the Gaussian filter, / w (t) = 

j+infty ^( T ) e -i™-c 2 (r-t) 2 A 2 dr , Daghed 

curves: the odd harmonics frequency multiplied by the factor (1 + /3)/(l — 0) 
calculated from the fast ion spectrum maximum at the time of reflection. Modes aliasing occurs at later times due to the fixed 
width of the filter and a fast change of the frequency multiplication factor. 
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has size of 50 A with the resolution of 128 steps per A along the a;-axis and 32 A with the resolution 16 steps per A 
along the y-axis. The number of quasi-particles is 10 6 . We note that the use of a circularly polarized driver pulse 
may provide a smoother start of the slab acceleration in the radiation pressure dominant regime [TT] , nevertheless it 
was chosen to be p-polarized in order to easily distiguish between the driver and the source pulses. In addition, our 
choice demonstrates the robustness of the double-sided mirror effect. The results are shown in Figs. [2]3] where the 
spatial coordinates and time units are in the laser wavelengths and wave periods, respectively. 

The driver laser makes a cocoon where it stays confined, Fig. |2ja). At t — 37 x 2n /uj, the ions are accelerated 
up to 2.4 GeV while the majority of accelerated ions carry the energy about 1.5 GeV, Fig. [2jb) . The accelerating 
plasma reflects the source pulse, which becomes chirped and compressed about 10 times, Fig. [2ja). As the mirror 
velocity, cj3, increases, the reflected light frequency grows as (1 + /3)/(l — (3), thus the electric field profile along the 
x-axis becomes more and more jagged, Fig. [3ja). A portion of the source pulse reflected from the curved edges of the 
expanding cocoon acquires an inhomogeneous frequency upshift determined by the angle of the reflecting region. At 
the begining, the magnitude of the reflected radiation is higher than that of the incident source (~ 3 times), due to the 
double Doppler effect and an enhancement of reflectivity owing to the plasma slab compression under the radiation 
pressure exerted by the driver pulse. In an instantaneous proper frame of the accelerating mirror, the frequency of 
the source pulse increases with time, thus the mirror becomes more transparent. Correspondingly, the source starts 
to be transmitted through the plasma more efficiently, as seen in Fig. [2ja) where the transmitted radiation is focused 
owing to a cocoon-like spatial distribution of the plasma. 

The reflected radiation has a complex structure of the spectrum, Fig. |3^b) . It contains not only the frequency- 
multiplied fundamental mode of the source pulse, but also high harmonics due to the nonlinear interaction of the source 
with the plasma slab. This is also seen in Fig. [2jc) , where the first two consecutive cycles of the reflected radiation 
exhibit presence of high harmonics, while the later cycle is compressed together with its harmonics in comparison with 
the earlier cycle. The spectrum is also enriched by a continuous component since the mirror moves with acceleration. 
The high harmonics generation efficiency is optimal for a certain areal density of the foil, according to the condition 
a s sa irn e lr e A s 7 , where r e = e 2 /m e c 2 is the classical electron radius. Initially far from this condition, the accelerated 
plasma slab satisfies it at certain time, when harmonics are generated most efficiently. 

In order to analytically describe the reflected EM wave we use the approximation of an infinitely thin foil (see 
also Refs.[H [7]), representing a mirror moving along the cc-axis with the coordinate Xu(fy- We consider the one- 
dimensional (ID) Maxwell equation 

d 2 *f c2 d 2 ^ | A7re 2 nJ5[x-^ M (t)} ^ = Q 
dt 2 ' dx 2 m-elM 

where jm — [l — {d3>M/dt) 2 c~ 21 \ is the Lorentz factor of the mirror, is the EM wave vectror-potential, 5 is 
the Dirac delta function. Let k is the incident wave number. Transfromations to dimensionless variables and to new 
variables, £, 77, which are the characteristics of the Maxwell equation, 

x = kx , i = kct , (2) 
£ = (i-*)/2, » 7 = (s + t)/2 J (3) 

jr M (t) = k-'x M {n - , *f{x, t) = ^A(t ri) , (4) 

and the property 8{kz) = k~ l 5{z) yield the equation 
where x = ^n e lr e \, A = 2ir/k, and 

4>fori) = t + v-XM(v-S) , (6) 

l^,r } )=[l-X' 2 {v-0V 2 ■ (7) 
We seek the solution to Eq. (|5| in the form of the incident, transmitted and reflected waves: 

^•"-{^), + ^)<^"" >0! < 8 > 

Here the factor e 2lri — e lk ( x + ct ) represents the incident wave. The solution should satisfy the boundary conditions 
at the position of the mirror, ij){£,ri) — 0. We introduce new functions £o(rf) and f]o(Oi which satisfy the following 
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expressions 



^(£0(77), 7?) = for V77, 
V>(f,rio(f)) = OforVf. 



(9) 
(10) 



The requirement that the solution is continuous, ^4(f, 770(f) — 0) = A(f, 779(f) + 0), leads to the following condition: 

ai (0 + a e 2 ^ =a 2 ( Vo (0). (11) 

The remaining conditions can be obtained from Eqs. ((5J and ((8). Integrating Eq. ^ over 77 in the vicinity of 770(f) 
for fixed f and some small e > 0, we obtain: 



dA 
5f 



(12) 



If (770(f) — £ ) 7^ 1 for given f , which means that the mirror velocity does not reach the speed of light in vacuum, 
the function ip has a simple zero at (f , 770(f)), therefore the equation 



1,0+6 f(f 770) 

<m r?)*7- - ' 



dtp/dr] 



(13) 



where the derivative dtp/drj is taken at the point (f, 770(f)), holds for any integrable function /. Using Eq. (13l we 
find that at the limit e — > Eq. ( 12 ) gives the magnitude of the jump discontinuity of the derivative A^ = dA/d£ at 
V = Vo(0 for nxed f : 



dA 

where we introduce the digamma factor, F , 



»7=no(«)+o 



xr(f,7? (f))A(f, 770(f)) 



»7=77o(S)-0 



1 + X' m (ti-0 
1 - X' M (77 - f ) 



1/2 



(14) 



(15) 



Similar expression is obtained for the magnitude of the jump discontinuity of the derivative A v — dAjdr\ at f = foC 7 ?) 
for fixed 77. For the ansatz ^ these expressions give the following two ordinary differential equations (ODE): 



a'i(f) = x(oi(0 + a e 2 ^) !*,(£)) , 



2ia e 2171 -a' 2 (r)) = 



02(17) 



(16) 
(17) 



The reflected, ai(f), and the transmitted, 0,2(11), waves are determined by Eqs. (Ill, (16), and ( IT I which can be 
easily reduced to quadratures. 

In the simplest case of uniform motion, 



^m(^) — & — const, 



have 



F(Z,v) = F = 



1 + /31 1/2 
1-/3. 
%(f) = 



The solution to Eqs. (Ill, (16 1, and (17 1 reads 

ffl! = - 



X + 27F 

«2 = ao — expiry), 
X + 27F 



(18) 

(19) 
(20) 

(21) 
(22) 
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so that the reflection coefficient in terms of the number of photons is 

2 / , M 2 



R 



(n e lr e \y 



{n e lr e \) 2 +A 1 2 M 



(23) 



thus we recover the result of Ref . [I] . 

In the case of a mirror moving with a uniform acceleration gkc 2 , for simplicity we consider the particular trajectory 



2il/2 



Then we obtain 



XM&^g-'li + igt) 2 } 



F(ZMO) = (200 -1 , F(£ (v),v) = 2 OT , 



and the solution to Eqs. (11), ( 16 1, ( 17 1: 



ai (0-^(^) 2s r ^(2 l5 2 £)-\0 
aa(»?) = ^(-2i»?r*r 



a Q e 



2f';/ 



where T(a, zx, z "i) — J z ^ t a * e t dt is the generalized incomplete gamma function |13j . At £ — > 0, 



23 



2ff 2 ? 



(e + o(e 2 )), 



(24) 



(25) 
(26) 



(27) 
(28) 

(29) 



where T(z) is the Euler gamma function |13j . The frequency of the reflected radiation increases as as m the 
case of a perfect mirror of Ref. |14j . However, in our case the mirror reflectivity decreases with time. An observer at 
infinity (which corresponds to £ = 0) see the radiation with the frequency spectrum of the following intensity 



^(l + §;2,2 + §;-^ 



(30) 



v 2 5 4(l + 2g/ X ) 2 

where v is the observed frequency, iF 2 (ai; bi, 62; z ) is the generalized hypergcometric function 13J. Here l v is defined 
as the square of the modulus of the Fourier transform of the function ai(£), I v = -7= e w ^ai(£)d{; 

cast out the essential singularity at v — representing the finite limit ai(£ — ► +00) 
intensity decreases with frequency, v, as 



where we 
For large v, the spectral 



4 5 4 



r(-£) 



X 
2.9 



cos 



(^) + sin(^) 



2ff 2 



2\ 5/4" 



(31) 



Now we consider the case of a mirror oscillating with frequency (normalized on the incident wave frequency) 



d ( m 

choosing the following trajectory of the mirror 



= <?cos(f2f) , 



„ 1 / cos(fii) 
Am i = — arctan , 

where /i 2 = 1 + il 2 /g 2 . Solving the equation </>(£, 770 (0) = 0, we obtain 



Vo(0 = ~ ^ arctan H - ~ 



fi£ + arctan h 1 

tT~ 2 



/itan(fif) + 1 



(32) 
(33) 

(34) 
(35) 
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where the function \_z] gives the integer closest to z. The digamma factor, Eq. (15 1, for 77 = ?7o(£) reads 



h 2 - 1 



-,1/2 



h 2 + l + 2hsm(2fl£_) 



(36) 



For Eq. (33 1 the only bounded solution to Eq. (16 1 is 




(37) 



£ W= ex p{ 9 (fe) F (^-fl(w)} 



(38) 



where F(z|m) is the elliptic integral of the first kind with an asymptotic oc z for z — > 00 |13j . 

In conclusion, a solid density plasma slab, accelerated in the radiation pressure dominant regime, can efficiently 
reflect a counter-propagating relativistically strong laser pulse (source). The reflected electromagnetic radiation 
consists of the reflected fundamental mode and high harmonics, all multiplied by the factor w 4j 2 , where 7 increases 
with time. In general, the reflected radiation is chirped due to the mirror acceleration. With a sufficiently short source 
pulse being sent with an appropriate delay to the accelerating mirror, one can obtain a high-intense ultra-short pulse 
of X-rays. 

For the mirror velocities greater than some threshold, the distance between electrons in the plasma slab in the 
proper reference frame becomes longer than the incident wavelength. Thus the plasma slab will not be able to afford 
the reflection in a coherent manner, where the reflected radiation power is proportional to the square of the number 
of particles in the mirror. In this case the reflected radiation becomes linearly proportional to the number of particles. 
Even with this scaling one can build an ultra-high power source of short gamma-ray pulses, when the interaction of 
the source pulse with a solid-density plasma is in the regime of the backward Thomson scattering. 

We can estimate the reflected radiation brightness in two limiting cases. For 27 < (n e \ s s ) l / Q , the reflection is coherent 
and for the reflected photon energy tuv the brightness is Bm ~ £ s (fiw) 3 A s /47r 5 7i 4 c 3 , where £ s is the source pulse energy. 
For larger 7, the interaction may become incoherent and should be described as Thomson scattering, which gives Bj* ~ 
a d £ s (fia;) 2 r e A2/87r 4 ft 3 c 2 A^ For £ a = 10J, A s = 0.8/Ltm, hw = IkeV (7 = 13), B M = 0.8 x 10 40 photons/mm 2 mrad 2 s, 
which is orders of magnitude greater than any existing or proposed source, |15j . For the same parameters and 
A rf = 0.8^m, a d = 300, hw = lOkeV (7 = 40), B T = 3 x 10 32 photons/mm 2 mrad 2 s. 

Employing the concept of the accelerating double-sided mirror, one can develop a relatively compact and tunable 
ultra-bright high-power X-ray or gamma-ray source, which will considerably expand the range of applications of the 
present-day powerful sources and will create new applications and research fields. Implementation in the "water 
window" will allow performing a single shot high contrast imaging of biological objects [16 . In atomic physics and 
spectroscopy, it will allow performing the multi-photon ionization and producing high-Z hollow atoms. In material 
sciences, it will reveal novel properties of matter exposed to the high power X-rays and gamma-rays. In nuclear physics 
it will allow studying states of high-Z nucleus. The sources of high-power coherent X-ray and ultra-bright gamma-ray 
radiation also pave the way towards inducing and probing the nonlinear quantum electrodynamics processes. 
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